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Abstract 

Astrophysical systems differ often in two points from 
classical thermodynamical systems: 1.) They are open and 
2.) gravity is a dominant factor. Both modifies the homo- 
geneous equilibrium structure, known from classical ther- 
modynamics. In order to study the consequence for struc- 
ture formation in astrophysical systems, we carry out N- 
body simulations of self-gravitating systems, subjected to 
an energy-flow. The simulations show that physically re- 
alistic, time-dependent boundary conditions can maintain 
a molecular cloud in a statistically steady state, out of 
thermodynamic equilibrium. 

Moreover we perform some simple "gravo-thermal" N- 
body experiments and compare them with theoretical re- 
sults. We fin d negative specific heat i n an energy range 
predicted by Follana & Lalicna (1999). 
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1. Introduction 

An energy-flow through a system is related to an entropy- 
flow. If the entropy-flow leaving the system is larger than 
those entering, then the system evacuates it's internal, by 
irreversible processes produced entropy to the outer world. 
Consequently order is created inside the system, which 
diverges from a classical thermodynamical equilibrium in 
a closed system. Far from equilibrium, the system is no 
longer characterized by an extremum principle, thus los- 
ing it's stability. Therefore perturbations can lead to long 
range order, through which the system acts as a whole. 
Such a behavior is well known in laboratory hydrodynam- 
ics and chemistry. The underlying concepts such as "dissi- 
pative structures" and 11 self- organization" were extensively 
studied (see e.g. Nicolis & Prigogine 1977). But despite 



of their popularity they are until now only little stud- 
ied in the context of non-equilibrium structures in self- 
gravitating astrophysical systems. Therefore we take up 
some ideas of these concepts and build a simple model of 
an open self- gravitating system. With this model we want 
to check if an energy-flow can maintain a self-gravitating 
system in an statistically stable state, out of thermody- 
namical equilibrium and if gravitation in combination with 



an energy flow can create structures with a higher degree 
of order. 

The simplicity of the model permits, furthermore, to 
perform some "gravo-thermal" experiments whose results 
can be compared with theory. 

2. Why Dissipative N-Body Systems? 

The use of hydrodynamics to simulate self-gravitating, 
clumpy gas may be problematic, due to the following rea- 
sons: 

1. ) During recent years observations revealed that molec- 
ular clouds are structured down to the smallest resolvable 
scales. These structures appear fractal. As a consequence 
density, temperature and related fields are not everywhere 
differentiable, rendering the use of hydrodynamics prob- 
lematic. 

2. ) Hydrodynamics is based on the assumption of a local 
thermodynamic equilibrium (LTE). But such an equilib- 
rium may not be established in gravitationally unstable 
media because the speed of matter disturbances is com- 
parable to the sound speed ( Pfenniger 1998 ). 

3. ) Thermodynamical formalism incorporating gravity yield 
results differing from those of classical thermodynamics. 
An example of this is the so called gravo-thermal catastro- 
phe. The application of these formalism lead among other 
things to spatially inhomogeneous equilibrium states. In 
hydrodynamics gravity is inserted in the Euler equation as 
an external force. But the thermodynamical variables are 
determined by a theory not taking gravity into account. 
Because gravity can change thermodynamical behavior, 
it is not a priori clear that hydrodynamical methods can 
simulate self-gravitating gas correctly. 

Bearing in mind these considerations we try an other 
approach in order to model self-gravitating interstellar 
gas. That is, we use dissipative particles, representing dense 



cloud fragments, to simulate cosmic gas (Huber & Pfcn- 
aiger 200 la|). 



Furthermore, with such a realization we can check some 
thermodynamic results of self-gravitating systems with 
softened potentials. 

3. Model 
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Figure 1. The potential as a function of r in units of 
Gm/R, where m is the particle mass. For £ > 1/3 the 
potential is repulsive on the softening length scale. 




Figure 2. Perturbation potentials as a function of time and 
ID-space. The space coordinate r is determined through 
r 2 = x 2 + y 2 + z 2 and x — y = z. The time is in- 
dicated in units of dynamical times. We show two per- 
turbation potentials with different frequency and ampli- 
tude. Full resolution version of this figure is available at: 
http://o bswww. unige.ch/ Preprints / dyn_ art. htm< 



where R is the distance from the center of the system. 



3.1. Principle 

To prevent gravitationally unbound particles from disso- 
lution, the particles are confined in a spherical potential 
well. This prevents a matter flow. However, our system is 
subjected to an energy flow. This flow is maintained by 
energy injection (heating) and dissipation. The energy is 
injected due to time and space dependent potential per- 
turbations. The dissipation is realized by a local or global 
cooling scheme. 

Among others we apply a particle potential becoming 
repulsive on the softening length scale. With such a po- 
tential we want to prevent the formation of very dense 
particle agglomerates below the resolution scale. The re- 
pulsive force is analogous to the molecular van der Waals 
force. 



3.2. Self- Gravity and Repulsion 

The gravitational forces are computed on the Gravitor 
Beowulf Cluster^] with a parallel tree code, based on the 
Barnes fc Hut (1986) tree algorithm. The particle poten- 
tial is: 
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(1) 



where m is the particle mass, e is the softening length and 
r is the distance from the particle center. The parameter 
£ determines the deviation from the Plummer potential 
$pi. The potential is repulsive for £ > 1/3 in the range 
M < eVT 3 ^ - !) ( see Figure |). 

3.3. Confinement 

The confinement potential, preventing the dissolution of 
gravitationally unbound particles has the form, 

gconf (X R 16 , (2) 

1 http:/ /obswww.unige.ch/~pfennige/gravitor/gravitor.html 



3.4. Perturbation (Heating Scheme) 

The energy injection is due to time-dependent boundary 
conditions or more precise due to a perturbation poten- 
tial. The perturbation potential is shown in Figure [| The 
amplitude and the frequency can be controlled by param- 
eters. If our system represents a molecular cloud then po- 
tential perturbations can be due to star clusters, clouds 
or other high mass objects passing in the vicinity. Indeed 
such stochastic encounters must be quite frequent in galac- 
tic disks and we assume, 
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(3) 



where / per t is the mean frequency of the encounters. Thus 
these encounters can provide a continuous low frequency 
energy injection on large scales. 

3.5. Energy Dissipation 

In order to maintain an energy-flow we must dissipate 
somehow the energy injected at large scales. We use two 
different dissipation schemes. 

1. Local Dissipation 

A particle dissipates energy during an "inelastic encounter" 
with an other particle. Thus we add friction forces, de- 
pending on the relative particle velocities V and positions 
r. The friction forces are: 



Fi cx 



V cos # , 







r > Ae 
r < Ae 



(4) 



where Vrcosfi = Vr, i = {x,y, z}, e is the unity vector 
V /V and A is a parameter > 1 ensuring the local nature 
of the dissipation. 
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Figure 3. A self- gravitating system subjected to an energy-flow. The energy-injection compensates the dissipated energy 
and a statistically steady state is reached. 
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Figure^. The evolution of a self- gravitating system with local dissipation, but without heating, i.e., without energy-flow. 



2. Global Dissipation 

The global dissipation term depends not on the relative 
velocity V but on the absolute velocity v. The global fric- 
tion force, causing an energy dissipation is: Fi cx —V{. 



4. Results 



We present here some preliminary results of a paper in 



preparation (Huber & Pfennigcr 2001b) 



4.1. Self-Gravity and Energy Flow 

With an appropriate tuning of the energy injection we 
can prevent the collapse of a self-gravitating dissipative 
system and maintain it in a statistically steady state out 
of thermodynamic equilibrium. This is shown in Figure ^|. 
The system subjected to an energy flow attains a statis- 
tically steady state after about 10rd yn - The evolution of 
the same system, but without energy injection, is shown 
in Figure ^. Both simulations were carried out with a local 
dissipation scheme and 10000 particles, from which 5000 
are shown. 

The system subjected to an energy-flow is not ther- 
malized. A temperature gradient arises, in such a way that 
a cooler mass condensation is embedded in hotter shells 
(see Figure ||) . Such temperature gradients are typical for 
molecular clouds. 

We checked the model behavior in dependence of the 
different particle potentials, heating and cooling parame- 
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Figure 5. The evolution of the Lagrangian temperature of 
the system shown in Figure [J. (3 is the inverse dimen- 
sionless temperature (/3 = l/T). Each curve depicts the 
temperature of a mass fraction, contained in a sphere. 
The mass spheres are centered at the center of mass 
and the corresponding mass fractions are: AM/M — 
{5%, 10%, 20%,..., 80%, 90%, 98%}. Circles: AM/M = 
5%. Crosses: AM/M = 98%. 



ters. Actually we find statistically stable states, but these 
states are not endowed with higher degrees of order. 
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Figure 6. The Follona (n = 10) and the Plummer (e = 
0.05) potential in units of Gm 2 / R, where m is the particle 
mass. 

5. Gravothermal experiments 

If we switch off the heating process and use a global dis- 
sipation scheme we can cool the system and maintain it 
nearly in thermodynamic equilibrium. Thus we can per- 
form some simple "thermodynamic experiments" of self- 
gravitating N-body systems. Follana & Laliena examined 
the thermodynamics of self-gravitating systems with soft- 
ened potentials. They achieve a softening by truncating to 
n terms an expansion of the Newtonian potential in spher- 
ical Bessel functions (we will call this potential hereafter 
Follana potential). This regularization allows the calcula- 
tion of the thermodynamical quantities of self-gravitating 
systems. The form of their potential is similar to a Plum- 
mer potential with a corresponding softening length (see 
Figure ^|) . Thus we can check their theoretical results by 
using a Plummer softened potential. 

Because the collapsed phase is probably sensitive to 
the short distance form of the potential we expect quali- 
tatively a similar behavior only for the Plummer softened 
potential (£ = 0) and a more different behavior for the 
potentials with short range repulsion (£ > 1/3). Figure [7] 
shows the evolution of the inverse temperature as a func- 
tion of the total energy^] e for two simulations with £ = 
and £ = 2/3, respectively as well as the analytical curve 
calculated by Follana & Laliena. The figure confirms our 
expectations. We find for the Follana and the Plummer po- 
tential a negative specific heat in the same energy range. 
The system with a repulsive potential re-enters the zone of 
positive specific heat earlier than those without such a re- 
pulsion below the softening length. After entering the zone 
of negative specific heat a collapsing phase transition takes 
place, separating a high energy homogeneous phase from 
a collapsed phase. The collapsed phase resulting from the 
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Figure 7. Inverse temperature (3 versus energy e resulting 
from models with softened potentials. The solid line cor- 
responds to the Follana potential. The other curves depict 
the evolution of the simulated systems. Crosses: Repulsive 
potential (£ = 2/3). Dots: Plummer potential (£ = 0). The 
circle indicates the initial states of the simulations. 

N-body simulations shows no core-halo structure. Such a 
structure is only formed due to relaxation after the cooling 
has stopped. 

6. Conclusions 
For systems subjected to an energy flow we find: 

— Potential perturbations caused by astrophysical ob- 
jects passing in the vicinity of a self-gravitating sys- 
tem can compensate the energy loss due to dissipation, 
thus prevent the system from collapsing and maintain 
a statistical steady state. 

— The statistical steady state is out of thermodynamic 
equilibrium and consists of a dense cold core moving 
in a hotter halo. 

The results of the gravo-thermal experiments are: 

— The range of negative specific heat agrees for a Plum- 
mer softening with the theory of self-gravitating sys- 
tems with softened potentials. 

— A softened potential incorporating short range repul- 
sive forces reduces the range of negative specific heat. 
Systems with such repulsive forces re-enter the zone of 
positive heat capacity, thus becoming stable, at higher 
energies than those with "conventional" softening. 
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